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ABSTRACT 



The brightest continuum source in the Orion Molecular Cloud-3 region (OMC- 
3), MMS 6, was observed with the Very Large Array (VLA), the Nobeyama 
Millimeter Array (NMA), and the Submillimeter Array (SMA). Our data were 
supplemented by near- to mid-infrared archival data taken by Spitzer Space Tele- 
scope. The compact continuum source, MMS 6-main, was detected with an H 2 



- 2 - 



mass of 3.0 M with a size of 510 AU. Despite its compact and well condensed 
appearance, neither clear CO outflow, radio jet, nor infrared sources (at a wave- 
length shorter than 8 fim) were detected at MMS 6-main even with the present 
high-spatial resolution and high-sensitivity observations. The derived H 2 column 
density, 2.6xl0 25 cm -2 , corresponds to a visual extinction of ^^15000 mag., 
and the derived number density is at least two orders of magnitude higher than 
for the other OMC-2/3 continuum sources. The volume density profile of the 
source was estimated to have a power-law index of 2 or steeper down to a radius 
of ~450 AU. The time scale to form a protostar at the center or the time scale 
elapsed after its formation is estimated to be 830 to 7.6xl0 3 yr. This is much 
shorter than the typical lifetime of the Class 0/1 protostars, which is ~ 

10(4-5) 

yr, suggesting that MMS 6-main is probably in either the earliest stage of the 
proto-stellar core or in the latest stage of the pre-stellar phase. 

Subject headings: ISM: clouds — ISM: individual (OMC-2/3) — stars: formation 
— ISM: pre-/proto-stellar cores — ISM: molecules — radio lines: ISM 



1. INTRODUCTION 

Stars form in dense molecular cloud cores (n^lO 4-5 cm -3 ), which are the precursors of 
protostars. In order to understand the central steps of the star-forming process, it is essential 
to study the structures, kinematics, and physical conditions of the parental molecular core. 
These properties provide the initial conditions for star-formation. 

One crucial method to trace the evolutionary status of the parental molecular cores is to 
investigate their density structures. Theoretical studies by Shu et al. (1987) have proposed 
the collapsing core model known as the "inside-out collapse model". Initially hydrostatic 
cores, with the density distribution of p(r)ocr -2 of a singular isothermal sphere, will dynam- 
ically collapse from the center with an expansion wave which propagates outward at the 
sound speed, C s . This infalling envelope has the size of C s t at the given time, £, and shows 
a density structure of p(r)ocr -L5 . For the alternative semi- analytical isothermal similarity 
solution for the collapse (Larson 1969; Penston 1969), the compression wave propagates in- 
ward at the sound speed, reaching the center at t — 0. At t > 0, this wave is reflected as an 
expansion wave, and propagates outward with the sound speed, C s . The density structure of 
this collapsing core is the same as for the inside-out collapse case, although the absolute value 
of the density is higher than the equilibrium case (4.4 times denser for the Larson-Penston 
solution) . 
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A number of prestellar- and protostellar-cores have been studied using infrared ex- 
tinction maps and millimeter/submillimeter emission images toward low-mass star- forming 
regions (e.g., Ward-Thompson 1994; Ward-Thompson 1999; Andre et al. 1996; Bacmann 
et al. 2000; Shirley et al. 2000; Kandori et al. 2005; Motte et al. 2001; Looney et al. 
2000; Harvey et al. 2003a, b). It has been suggested that the observed density structures of 
low-mass protostellar cores are roughly consistent with those predicted by either the inside- 
out collapse model or the isothermal sphere model with the density profile p(r)^r -L5 or 
p(r)~r -2 . On the other hand, observed prestellar cores show a shallower density power-law 
index (i.e., flattened density profile) at the core center. This can be explained if (i) the core 
is not sufficiently evolved to have reached the isothermal sphere density profile (i.e., t<0 case 
in the Larson- Penston case), (ii) other mechanism such as a magnetic field could support the 
system (e.g., Crutcher et al. 1994), or (hi) the core has a density structure explained by a 
finite-sized Bonner-Ebert sphere (e.g., Bonnor 1956). A large sample of these studies is one of 
the crucial strategies for understanding the initial conditions of the star-forming core as well 
as investigating the evolutionary status of the central stars. In particular, high-resolution 
observational studies with a spatial resolution of much less than 1000 AU are necessary to 
constrain the radius where the density profile steepens and to compare this with the theo- 
retical model for the very early stages of protostellare cores. Such higher angular resolution 
studies are still limited in number. 

The OMC-2/3 region is located at the northern-part of the Orion Molecular Cloud 
A (rf-400 pc; Sandstrom et al. 2007; Menten et al. 2007; Hirota et al. 2007; Kim et 
al. 2008), which is the nearest giant molecular cloud, and one of the best-studied cluster- 
forming regions at all observable wavelengths (e.g., Bally et al. 1987, Allen et al 2007; 
Tatematsu et al. 1993, Johnstone & Bailey. 1999; Tsujimoto et al. 2005). In this region, 
filamentary-shaped molecular clouds, lying in the north-south direction with several pc scale, 
were discovered in the submillimeter continuum observations with the IRAM 30 m telescope 
and the JCMT/SCUBA (Chini et al. 1997; Johnstone & Bally. 1999; Johnstone et al. 2003). 
These filamentary structures contain a few tens of millimeter and sub-millimeter continuum 
sources, presumably tracing prestellar- and protostellar-cores (Chini et al. 1997; Lis et al. 
1998; Nielbock et al. 2003; Takahashi et al. 2006; Takahashi 2008a, 6; Shimajiri et al. 2008). 
In addition, approximately a dozen radio jets, molecular outflows, and shock excited H 2 
emissions, were detected in this region (Yu et al. 1997; Reipurth et al. 1999; Aso et al. 
2000; Stanke et al. 2002; Williams et al. 2003; Matthews et al. 2005). Moreover, our recent 
extensive molecular outflow survey with ASTE in the CO (3-2) emission have revealed 14 
molecular outflows, which includes seven newly detected outflows, suggesting that the OMC- 
2/3 region is one of the most active nearby on-going star- forming regions and contains the 
richest assemblages of protostars (Takahashi et al. 20086). 
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From a series of our OMC-2/3 survey observations combined with previous millimeter 
and submillimeter continuum studies, we have found an interesting 1.3 mm continuum source, 
MMS 6, which is the brightest source in the millimeter- and submilimeter- wavelengths in the 
OMC-2/3 region. This brightest source is located at the center of the OMC-3 region, and the 
1.3 mm flux is roughly one order of magnitude larger than those for any other continuum 
sources in OMC-2/3 (Chini et al. 1997; Johnstone & Bally 1999). Despite the unusual 
appearance of MMS 6, no signature of star- format ion activities such as molecular outflow 
or jet have been detected toward MMS 6 by Takahashi et al. (20086). The bolometric 
luminosity, temperature, and core mass of MMS 6, which were derived from the previous 
single-dish millimeter to submillimeter observations combined with IRAS data, are <6OL , 
15-25 K, and 36 M , respectively (Chini et al. 1997). The derived luminosity and mass 
are one order of magnitude higher than the values derived in typical low-mass star-forming 
regions such as the Taurus molecular cloud (e.g., Myers & Benson 1983). 

That MMS 6 does not show any significant signatures of star- format ion activities, sug- 
gests that MMS 6 is at the earliest stage of evolution. To address this issue, our primary 
goal is to constrain the physical conditions such as size, mass, density, density profile, and 
dust grain properties. Comparison of the derived values with the low-mass core values is 
also important to elucidate the star- formation environment. There is another close- by 1.3 
mm source, MMS 5, which is located at 25" north-west of MMS 6. We have detected a clear 
bipolar outflow associated with MMS 5 in the CO (3-2) emission, which is mainly elongated 
in the east-west direction. While there may be additional diffuse gas extending to the MMS 
6 direction, the relatively large beam of the previous CO (3-2) observations (26") makes it 
difficult to investigate the detailed spatial structures of molecular outflows and the physical 
properties of cores in the MMS 5/6 regions. 

To unveil the evolutionary status of MMS 6, in this paper, we have investigated the 
detailed physical properties with higher spatial resolution observations using the VLA, the 
NMA, and the SMA. Moreover, near- to mid-infrared data retrieved from the Spitzer Space 
Telescope (IRAC and MIPS) archival data, and near-infrared JHK S data taken with the 
SIRIUS/IRSF (taken by Takahashi et al. 20086), were compared with the interferometric 
data with relative position accuracy. 

The details of the millimeter- and submillimeter-interferometric observations are pre- 
sented in Section 2. Results of the millimeter- to infrared-continuum emissions and also 
some molecular lines are shown in Section 3. Physical properties of MMS 6 are discussed in 
Section 4. Finally, Section 5 summarizes the paper and shows future prospects. 
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2. OBSERVATIONS AND DATA REDUCTION 

2.1. Interferometric Observations 

Aperture synthesis observations of molecular lines and continuum emissions were carried 
out using the Nobeyama Millimeter Array (NMA), the Submillimeter Array (SMA), and the 
Very Large Array (VLA). The phase reference center was set on the peak position of the 1.3 
mm source, MMS 6, which has a coordinate of R.A. (J2000) = 5 /l 35 m 23.48 5 , Decl.(J2000) = 
-05°01 / 32 // .20 (Chini et al. 1997). The observed parameters for the continuum and molecular 
line observations are summarized in Table 1 and Table 2, respectively. 

2.1.1. NMA Observations 

The 3.3 mm and 2.3 mm continuum emission were observed with the six 10 m antennas 
of the Nobeyama Millimeter Array (NMA) 0. H 13 CO + ( J= 1-0; 86.754 GHz) emission was 
observed together with the 3.3 mm continuum emission, and 12 CO (J— 1-0; 115.271 GHz) 
emission was observed independently. These observations were performed during May 2004 
to January, 2006. The 3.3 mm and 2.3 mm continuum data were obtained using a digital 
wideband correlator, UWBC, which has a bandwidth of 1024 MHz and 512 frequency chan- 
nels for each of the lower- and upper-sidebands (LSB and USB; Okumura et al. 2000). The 
LSB and USB continuum data separated by 12 GHz in RF were combined to improve the 
sensitivity. Spectral information of these lines were obtained with the 32 MHz bandwidth 
FX correlator containing 1024 channels. The rest frequencies of the H 13 CO + and 12 CO lines 
were placed in the center of the LSB and USB, respectively. The phase and amplitude cal- 
ibrator, 0531+135, was observed every 20 minutes. Observations of Uranus provided the 
absolute scale for the flux density calibration. The overall flux uncertainty was estimated 
to be ~15%. The passband across the bandwidth was determined from observations of 3C 
454.3 at 3.3 mm and 0530+135 at 2.3 mm with a 30 minute integration. The combined 
arrays provided projected baselines ranging from 3 to 115 kA at 3.3 mm and 4 to 38 kA at 
2.3 mm. 

The raw data were calibrated for the complex gain and passband using the standard 
reduction program, UVPROC-II (Tsutsumi, Morita, & Umeyama 1997). After the calibra- 
tions, final CLEANed images were made using the Astronomical Image Processing System 



1 Based on the observations made at the Nobeyama Radio Observatory, which is a branch of the National 
Astronomical Observatory, an inter university research institute operated by the Ministry of Education, 
Culture, Sports, Science, and Technology of Japan 
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(AIPS), developed at NRAC0. Natural weighting and robust weighting (0 = the weighting 
between natural and uniform) were applied for the line data and continuum data, respec- 
tively. The resulting synthesized beam sizes, the achieved rms noise levels as well as other 
observational parameters are summarized in Table 1 and Table 2 for the NMA molecular 
lines and continuum observations , respectively. 



2.1.2. SMA Observations 

The 0.9 mm continuum data were taken using the eight 6 m antennas of the Submil- 
limeter Array | (SMA; Ho et al. 2004), on January 2 and October 21, 2008. The combined 
arrays provided projected baselines ranging from 13 to 88 kA. Both the LSB and USB data 
were obtained simultaneously with the 90° phase switching technique by the digital spectral 
correlater, which has a bandwidth of 2 GHz in each sideband. The LSB and USB continuum 
data separated by 10 GHz RF were combined to improve the sensitivity. The phase and 
amplitude calibrator, 0531+135 (1.7 Jy), was observed every 20 minutes. Observations of 
Titan provided the absolute scale for the flux density calibration. The overall flux uncer- 
tainty was estimated to be ~20%. The passband across the bandwidth was determined from 
observations of 3C 273 and 3C 454.3 with a 60 minute integration in each. 

The raw data were calibrated using MIR, originally developed for the Owens Valley 
Radio Observatory (Scoville et al. 1993) and adopted for the SMA. After the calibration, 
final CLEANed images were made using the AIPS task IMAGER with a robust weighting 
(=0). The resulting synthesized beam size was 2 // .2x2 // .0 (or 880 AU x 800 AU at the 
adopted distance of 400 pc) with a position angle of -62°. The achieved rms noise level was 
43 mJy. The observational parameters are summarized in Table 1. 



2.1.3. VLA Observations 

The 3.6 cm data were taken with the Very Large Array (VLA) of NRAO in the CnD- 
array on June 1, 2008. Both the LSB and USB data were obtained simultaneously with 



2 NRAO is a facility of the National Science Foundation, operated under cooperative agreement by Asso- 
ciated Universities, Inc. 

3 The Submillimeter Array is a joint project between the Smithsonian Astrophysical Observatory and the 
Academia Sinica Institute of Astronomy and Astrophysics and is funded by the Smithsonian Institution and 
the Academia Snica. 
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the 90° phase switching technique by the digital spectral correlater, which has a bandwidth 
of 50 MHz in each sideband (total bandwidth of 100 MHz). The array provided projected 
baselines ranging from 8.7 to 55 kA. The phase and amplitude calibrator, 05416-05418 (1.4 
Jy), was observed. The 7.3 mm data were also taken with the VLA in the CnD- array. 
In addition, C-array data retrieved from the NRAO archival system were combined. The 
combined arrays provided projected baselines ranging from approximately 4.5 to 477 kA. 
The phase and amplitude calibrator, 05416-05418 (0.7 Jy), was observed. In order to track 
the short-term amplitude and phase variations induced by the earth atmosphere, the fast 
switching technique was adopted for both 3.6 cm and 7.3 mm observations. The cycle time 
of the observing target, MMS 6, and the phase-amplitude calibrator, 05416-05418 which lies 
about 2 ° from MMS 6, was set at 4.3 minutes. For the 7.3 mm observations, the pointing 
accuracy of each antenna was checked every 45 minutes from scans of the quasar 0541-056. 
Absolute flux calibrations at 3.6 cm and 7.3 mm were performed by observations of 3C 286 
for the C-array data and 3C48 for the CnD-array data. The overall flux uncertainties of 3.6 
cm and 7.3 mm observations were estimated to be a few percents. 

The raw data were calibrated using the AIPS software package. After the calibrations, 
we made final images using the AIPS task IMAGR with robust weighting (=0). The resulting 
synthesized beam sizes at 3.6 cm and 7.3 mm were 10 // .0x3 // .0 (or 4000 AU x 1200 AU at 
an adopted distance of 400 pc) with a position angle of 67° and l".4x0".6 (or 560 AU x 240 
AU) with a position angle of -87°, respectively. The achieved rms noise levels at 3.6 cm and 
7.3 mm were 50-60 /xJy and 19 mJy, respectively. Further 3.6 cm and 7.3 mm observational 
parameters are summarized in Table 1. 

3. RESULTS 
3.1. Source Identification Based on Continuum Images 

3.1.1. Millimeter and Submillimeter Continuum Emission 

Figure 1 shows the four-band continuum maps (i.e., 0.9, 2.3, 3.3, and 7.3 mm) toward 
MMS 6 taken with the VLA, NMA, and SMA. The single continuum source detected by the 
single-dish 1.3 mm, 850 /xm, and 350 /xm observations (Figure la), MMS 6, was resolved to 
be a strong main component (hereafter MMS 6-main) with additional substructure by the 
interferometric observations. To measure the size, position, and flux, 2D-Gaussian fitting 
was applied to each map, and the derived parameters are summarized in Table 3. 

MMS 6-main, which has the peak position of [RA(J2000)=5 /l 35 m 23.43 s , Dec(J2000)=- 
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5°01 / 30.6 // ], was detected at all the observed millimeter- and submillimeter- wavelengths. The 
peak position of the 7.3, 3.3, 2.3, and 0.9 mm data coincides within the positional accuracies 
of these data, <0.5". These peak positions are located at approximately 2" north-west of 
the 1.3 mm continuum emission peak, measured with the IRAM 30 m telescope (Chini et al. 
1997), which has a positional accuracy of ±5". The flux density of MMS 6-main measured 
with the SMA 0.9 mm continuum observations is 4.0 Jy, corresponding to 60 % of the flux 
measured with the 13 " beam of the JCMT/SCUBA 850 fim observations, which is 6.6 Jy 
(Johnstone & Bailey. 1999). 

A second component, which is located at ~1.5" north-east of the MMS 6-main (hereafter 
MMS 6-NE), was clearly detected at 0.9 mm. The other three- wavelength data also show 
a hint of the additional component with 3 to 6 a level emission. 2D-Gaussian fitting was 
applied to this second component in the 0.9 mm data. The peak position at the 0.9 mm 
emission coincides with the previously detected 3.6 cm source position by Reipurth et al. 
(1999) with a positional offset of 1.0", which is comparable to the absolute positional accuracy 
of the 3.6 cm observations (~1")- 

Extended emission from the large-scale filamentary structure, which was detected in the 
millimeter- and submillimeter- single-dish observations (Figure la), was not detected in the 
interferometric observations. As described in Table 1, the maximum detectable angular sizes 
for the interferometric observations, were 13 // -56 // . Therefore, the structures, which have a 
size scale larger than the maximum detectable sizes were missed in our images. 

In addition, another 1.3 mm source, MMS 5 identified by IRAM 30 m observations 
(Chini et al. 1997) was detected at approximately 22" north-west of MMS 6-main with the 
7.3 mm and 3.3 mm continuum emissions with 3-6 a level signals (Figure lc and d). 

3.1.2. Near- to Mid- Infrared Sources 

Figure 2a, b show the 2.2 /xm and 8 /im images overlaid with the 3.3 mm continuum 
emission. From the images, a total of five infrared sources located on the main OMC filament 
detected at either 2.2 /xm-, or 8 /xm-bands, or both (i.e., denoted by Figure 2 doted circles) 
were identified within ^0.5 pc radius from MMS 6-main (see Table 4). From the SEDs and 
color-color-diagrams, three of the sources (IRS 1, IRS 3, and IRS 5) have been identified as 
protostars (Peterson & Megeath 2008; see also Table 4). 

Of the five detected infrared sources, IRS 3 is the only one closely associated with the 
strong continuum source, MMS 6. The position of IRS 3 was measured in the Spitzer/IRAC 
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images 0. MMS 6-main and MMS 6-NE, are located at 2".l and ;/ .8 away from IRS 3. 
The important question is whether IRS 3 is a heating source of the millimeter/submillimeter 
sources, MMS 6-main or MMS 6-NE, or not? A comparison of the relative positions between 
the infrared data and the interferometric data in the context of the accuracy of the positional 
alignment between the images is necessary. The positional alignment of the near- to mid- 
infrared images, taken by SIRIUS/IRSF and Spitzer/IRAC, is determined by the standard 
stars in the observed frames, and the positional accuracy of these data is determined by the 
positional alignment accuracy of the 2MASS point source catalog. In this case, both near- 
and mid-infrared images have the same positional accuracies with a typical value of < ±l 7/ .2. 
The relative positions of the interferometric data are determined by the positional accuracy 
of the gain calibrators, which typically is < ±0".5. The best way to align the infrared and 
interferometric images is to match the common sources in the observed field. However, no 
common standard star was identified within our observed fields of view. Instead, we use 
another nearby OMC-3 millimeter source, MMS 7, observed by Takahashi et al. (2006), 
which was detected in both the interferometric 3.3 mm data taken with the NMA and the 
8 fim data taken with Spitzer/IRAC. In the 8 fim data, both MMS 7 and MMS 5/6 are 
in the same processed image. The 3.3 mm observations of MMS 7 and MMS 5/6, while 
not in the same primary beam, were made during the same 2004 November to 2005 March 
period with same array configurations. Since the positions of MMS 7 at 3.3 mm and 8 fim 
coincide to within 0.31", the accuracy of the positional alignment between the infrared- and 
interferometric images is on this order. The positional offset between MMS 6-main and IRS 
3, at 2.1 // , is formally a 7a result. This suggests that the infrared source, IRS 3, is likely not 
associated with MMS 6-main. On the other hand, the peak position of MMS 6-NE which 
was detected in the 0.9 mm continuum image, coincides with IRS 3 within the positional 
accuracies. This suggests that IRS 3 may be associated with the dusty envelope detected in 
the 0.9 mm emission. 

Another infrared source, IRS 1, is located at the position of the 1.3 mm continuum 
source, MMS 5 which also has been detected in the interferometric observations at 3.3 mm 
and 7.3 mm. The infrared source was only detected in the wavelength longer than the 
Spitzer/IRAC 5.6 //m-band, and the source has a Class type SED. 



4 We found that the peak position of IRS 3 as measured in the Ks image is roughly 1.5" shifted to the west 
from the peak position as determined from Spitzer/IRAC band images, which all self consistent to within 
^0.5". We have checked the positional alignments of the other point-like source in the same frame between 
different data sets (i.e., comparison of the SIRIUS/IRSF and Spitzer/IRAC data based on the 2MASS data), 
and found the point-like sources in the same frame coincide with each other. Hence the IRS 3 positional 
offset between the Ks-band and Spitzer/IRAC images is considered to be real. The inward extinction of K s 
most likely does not allow the true position of IRS 3 to be detect at K s 
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3.2. Molecular Line Emissions 

3.2.1. Search for high-velocity gas components in the MMS 5/6 region; CO (1-0) Emission 

Figure 3 shows the blueshifted (Vlsr=8.3 - 10.3 km s _1 ) and redshifted (Vlsr— 12.5 
- 14.5 km s _1 ) ranges of the 12 CO(1-0) emissions towards the MMS 5 and MMS 6 region, 
superposed on the Xs-band image taken with the SIRIUS/IRSF by Takahashi et al. (20086). 

A pair of blueshifted and redshifted components are seen in the east-west direction at 
east of MMS 6-main with a velocity less than 3.5 km s _1 from the systemic velocity. 
Distribution of the redshifted emission is highly consistent with the jet-like feature detected in 
the X s -band image emanating from IRS 4. The extended blueshifted component is elongated 
along the northeast to southwest. The counterpart of the blueshifted component is not 
detected in the X s -band. 

The extended redshifted emission with a velocity of 14.3 km s _1 <Vlsr< 15.4 km s" 1 was 
also detected in the CO (3-2) emission by the single-dish ASTE 10 m telescope (Takahashi 
et al. 20086) with the same position as the CO (1-0) emission detected with the NMA. 
At the position of MMS 6-main, no clear signature of typical molecular outflows, which 
has highly collimated bipolar structure with high- velocity component, was detected in the 
interferometric CO (1-0) observations with these high-spatial resolution (i.e., ^3000 AU) and 
high-sensitivity (i.e., ~ 6.5xl0~ 4 M ) observations. 

Another red-shifted component was detected to the east of MMS 5. The redshifted 
emission in the velocity range of Vlsr = 13.2- 21.9 km -1 was already detected in the CO(3- 
2) emission with the ASTE 10 m telescope (Takahashi et al. 20086). The distribution 
of the redshifted component detected by the NMA is consistent with that found by the 
ASTE CO(3-2) observations. The blue-shifted component associated with MMS 5, which 
was detected in the ASTE CO (3-2) observations, was not detected in the CO (1-0) emission, 
probably due to the primary beam attenuation of the NMA. 

3.2.2. Dense Gas Distribution; H 13 CO + (1-0) Emission 

In Figure 4, H 13 CO + shows an elongated structure from northwest to southeast with a 
ridge-like structure from the MMS 5 to MMS 6 region. This elongation agrees well with the 
elongation of the large-scale filamentary structure detected by single-dish observations at 1.3 
mm, 350 /im, and 850 fim wavelengths (Chini et al, 1997; Lis et al. 1998; Johnstone & Bailey. 
1999). Assuming that the H 13 CO + emission is optically thin and that the distribution of 
the excitation temperature is uniform, the total LTE mass of the entire H 13 CO + emission 
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structure is estimated to be 62-153 M with an average column density of (0.7-1.0) xlO 24 
cm -2 . Here, the excitation temperature of 20 K (e.g., Cesaroni et al. 1994, Chini et al. 
1997) and the H 13 CO+ abundance of (0.45-1.4) xl(T 10 (Aso et al. 2000, Takahashi et al. 
2006) are adopted. 

The position of H 13 CO + emission associated with MMS 6-main is established with a 2D- 
Gaussian fit with a deconvolved size of ll // x7 // .l (P.A. 9.7°), which corresponds to 4400x2800 
AU at rf=400 pc. The LTE mass derived from this emission is estimated to be 7.9-25 M 
with an average hydrogen column density of (0.9-2.7) xlO 24 cm -2 , under the assumptions 
described in the previous paragraph. The peak position of the H 13 CO + emission has a slight 
offset (Aoffset=2 // .5) from the dust continuum peak (i.e., MMS 6-main). Although the 
physical conditions cannot be tightly constrained with only the H 13 CO + (1-0) observations, 
the peak position offset between the H 13 CO + emission and the continuum emission could 
be explained by an opacity effect, density differences along the line of sight, or H 13 CO + 
depletion in the densest region , as traced by the millimeter- and submillimeter- continuum 
emissions. 

The H 13 CO + emission was also detected at MMS 5 with a deconvolved size of IQ" .OxQ" .2 
(P.A. 180°), which corresponds to 6400x2500 AU at rf=400 pc. The LTE mass derived 
from the emission was estimated to be 7.4-23 M with an average column density of (0.7- 
2.2) xlO 24 cm -2 under the assumptions described in the first paragraph. The H 13 CO + peak 
coincides with an infrared source IRS 1 and also dust continuum emission detected by our 
interferometric observations, suggesting that IRS 1 (MMS 5) is still surrounded by the thou- 
sands AU scale dense gas envelope. No strongly concentrated H 13 CO + emission was detected 
at the other three infrared positions of IRS 2, 4, and 5. The extended H 13 CO + emission 
with approximately 15 a level emission was detected at IRS 3, although it does not coincide 
with any of the H 13 CO + emission peaks. IRS 3 probably is still surrounded by the dense 
envelope, which was detected in the H 13 CO + emission. 



3.3. Spectral Energy Distribution 

Figure 5 shows the flux density of MMS 6-main as a function of the frequency. In order 
to enhance the sensitivity to the compact structure associated with MMS 6-main, the flux 
densities were measured in images where the uv-iange is restricted to be lager than 30 kA. 
We used the natural weighting map to be more sensitive to the extended structure at 7.3 mm 
since the 7.3 mm continuum image shows a smaller size-scale structure. Flux density and 
deconvolved size of the MMS 6-main compact component at each wavelength are summarized 
in Table 5. 
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3.6 cm (i.e., 8.4 GHz) continuum observations were made by Reipurth et al. (1999) 
with the VLA-D array toward the MMS 6 region. They have reported the detection of the 
3.6 cm continuum emission with a flux density of 0.15 mJy. However, the peak position of 
the 3.6 cm continuum emission is 2.0" offset from the the peak position of MMS 6-main. 
This positional offset is larger than the positional accuracy of the 3.6 cm observations (^l")? 
suggesting that the 3.6 cm continuum emission is probably not associated with MMS 6-main. 
In addition, no 3.6 cm continuum emission was detected in our VLA-DnC array observations 
with the upper limit (3 a level) of 0.09 mJy. Because of these uncertainties, in this paper, 
the 3.6 cm flux density reported by Reipurth et al. (1999) is adopted as an upper limit. 

The best fit of the spectral index a (i.e., F„(xis a ) in Figure 5 is derived to be 2.93. 
Since the spectral index, a ) and the dust emissivity, /?, have a relation of a = /3+2, the dust 
emissivity, /?, is estimated to be 0.93. 3.6 cm and 7.3 mm continuum emission sometimes 
trace an ionized jet (i.e., free- free emission) driven by low- to intermediate-mass protostars 
(e.g., Anglada 1995; Rodriguez 1997). Hence, if excess emission at centimeter wavelengths 
(i.e., mainly A>7 mm for the low- to inter-mediate mass case) were detected, it may indicate 
the presence of a central heating source. However, no clear excess emission caused by the 
free- free emission is detected toward MMS 6-main (Figure 5). 

On the assumption that the dust emission is optically thin at 0.9 mm and the tem- 
perature distribution of the dust continuum emission is uniform, the mass of the dusty 
condensation is estimated to be 3.0 M with a size of 400x320 AU (see flux and decon- 
volved size at 0.9 mm in Table 5), and the molecular hydrogen column density is estimated 
to be 2.6xl0 25 cm" 2 . Here, we assume that /3=0.93, T d =20 K, and ^ A =0.037 cmV 1 (A/400 
(Cesaroni et al. 1994; Chini et al. 1997; Keene et al. 1982). The hydrogen col- 
umn density and the visual extinction have a relation of A v = 0.56N H + 0.23, where N H is 
the equivalent column density in units of 10 21 cm -2 (Predehl & Schmitt 1995). The visual 
extinction is estimated from the hydrogen column density to be approximately 15000 mag. 
Moreover, this visual extinction corresponds to A 8/JjTn =790 mag., which is much higher than 
the value toward YSOs in the OMC filament, from several to 12 mag. (Peterson & Megath 
et al. 2008). Therefore, it is most likely that the central protostar associated with MMS 
6-main cannot be detected because of the extremely high-extinction toward MMS 6-main. 

Similarly, H 2 mass of MMS 5 were derived as 0.3 M from 3.3 mm continuum flux with 
assumptions that /3=0.8, 7^=20 K, and ft A =0.037 cm 2 g" 1 (A/400 /im) _/3 . Here, (3 were 
derived using 3.3 mm and 7.3 mm continuum data. The visual extinction was estimated to 
be ~220 mag., which corresponds to Ag /2rn =12 mag. This value is comparable to detected 
YSOs in the OMC filaments (Peterson & Megath et al. 2008). The central infrared source, 
IRS 1, was actually detected in the 24 /xm Spitzer/MIPS image. 
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4. Discussion: Nature of Dusty Core: MMS 6-main 

We are studying the physical properties of MMS 6-main in order to address its evolu- 
tionary status, as an example of the star- forming environment of the OMC-2/3 region. In 
this section, we discuss the physical properties of MMS 6-main from three points of view; (i) 
dust grain properties, (ii) comparison of physical properties between MMS 6-main and the 
other pre- and proto-stellar cores in the OMC-2/3 region, and (hi) internal structure of the 
cores. 



4.1. Dust Grain Properties 

The dust emissivity spectral index, /?, is related to dust properties such as grain sizes 
and grain shapes (e.g., Ossenkopf & Henning 1994; Pollack et al. 1994; Draine 2006). 
Observations have shown that (3= 1.5-2 in the general interstellar medium or dense molecular 
clouds/cores (e.g., Goldsmith et al. 1997), and in circumstellar disks (e.g., Beckwith & 
Sargent 1991; Andrews k Williams 2007). 

The anomalously low dust emissivity index, ^=0.3, in the MMS 6 region were reported 
from previous JCMT/SCUBA observations by Johnstone & Bally (1999). A similarly low 
value, /3=0.2, was also derived in the Orion KL region, which is the most active star-forming 
region in the OMC filament. Johnstone & Bally (1999) found that only these two regions 
show such low-/? values, while the other eighteen bright sources in the OMC filament show 
an average of (3=1.2. In Section 3.3, the dust emissivity spectral index of MMS 6-main was 
estimated to be (3 =0.93 from the multi- wavelength interferometric data. The estimated f3 is 
higher than the value derived from low-spatial resolution JCMT/SCUBA data toward MMS 
6, suggesting that there is a distribution in (3 for different size-scales. On the other hand, 
the derived (3 shows a lower value as compared with the (3 derived from the other bright 
continuum sources in Orion. 

In order to consider the behavior of (3 , the spatial distribution of the dust opacity 
spectral index, /?, is calculated using the 3.3 mm and 0.9 mm continuum images with the 
following equation, 



where F v is the flux in each pixel, and v\ and isq are the frequencies corresponding to 0.9 
mm and 3.3 mm, respectively. For the calculation, only regions which have a signal-to-noise 
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ratio above 3 a were used. Figure 6 presents the spatial distribution of f3 in the MMS 6 
region. In Figure 6, there is a suggestion that the (3 value has significantly changed along the 
north-south direction with a range of /3=0-1.5. Especially, (3 =0-1.5 is suggested at MMS 
6-main, while /?~1.5 is suggested at MMS 6-NE. These results imply that MMS 6-NE and 
MMS 6-main could have different dust grain properties. In the map, unreasonably small (3 
values, (3 <0, are shown at the southern part of MMS 6-main, which may be due to the 
contamination from an extended component. Note that there is extended emission in the 
3.3 mm continuum image, which presumably comes from ambient gas along the same line of 
sight, but no extended emission was detected at 0.9 mm (see Figure 1). 

From the multi- wavelength data, it is suggested that both MMS 6-main and MMS 6-NE 
are still deeply embedded in the molecular cloud. Nevertheless, the derived value for (3 is 
low. Draine et al. (2006) shows that the (3 value is highly dependent on grain size; i.e., large 
grain size gives a smaller (3 value. The observed low-/? could indicate faster grain growth in 
the densest region of the molecular core. Actually, recent interferometric observations have 
also reported /3^0.5-l toward Class type stars such as L1448 IRS2, L1448 IRS 3, and 
LI 157 (Kwon et al. 2008), and those results also suggest the faster grain growth. 

In addition to the grain growth, an unresolved optically thick component would also 
affect the (3 value. The average flux density was derived as 3.4 Jy at 0.9 mm. This corre- 
sponds to an averaged brightness temperature of 42 K. Similarly, brightness temperatures 
were derived from other three-bands data as 12.7 to 15.4 K. Since the opacity, the dust 
temperature and the observed brightness temperature have a relation of T5 = T^(l — e~ Tv \ 
if dust temperature is assumed at the typical temperature of the OMC filaments, T^~20 K 
(e.g., Cesaroni & Wilson 1994), the optical depth of the continuum structure is estimated to 
be 0.6-2.1, which is partially optically thick. Hence, one possible interpretation to explain 
the low-/? is that the larger outer envelope is likely more optically thin and have a higher 
value of /?, while the inner unresolved component may be more optically thick resulting 
in a smaller estimated (3 value. Moreover, higher opacity, due to condensed asymmetrical 
structures such as circumstellar disks inside the cores, can also produce a smaller (3 than 
the actual (3 as determined by the dust grain properties (e.g., Nakamoto & Kikuchi 1999; 
Nakazato, Nakamoto, & Umemura 2003). 



4.2. Density and Gravitational Stability 

In the previous studies, several other compact continuum sources (i.e., 500-2800 AU) 
were also detected in the 3.3 mm continuum emission toward other OMC-2/3 sources with 
the NMA (Takahashi et al. 2006, Takahashi et al. 2008a, Shimajiri et al. 2008). In 
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order to compare physical properties of these dust continuum sources with MMS 6-main, 
the average molecular hydrogen volume density was estimated using the equation of n = 
M H2 / {^t{R z iittih) with the assumption of a spherically symmetrical structure. Here, M H2 
is the mass derived from the dust continuum emission, \i and ttih are the mean molecular 
weight of 2.33 and hydrogen mass, respectively. Adopting M H2 =3.0 M Q and i?~510 AU, 
which were derived in Section 3.2, the mean volume density of MMS 6-main is estimated to 
be 6.5xl0 9 cm -3 . 

Takahashi et al. (2006) have investigated another protostellar core, MMS 7. The 3.3 
mm continuum emission was detected in the MMS 7 region with a mass and radius of 0.36- 
0.72 M Q and i?^500 AU, respectively, and hence the average volume density of H 2 was 
estimated to be (2.6-5.2) xlO 7 cm -3 . In addition to MMS 7, the typical hydrogen volume 
density toward the other dust continuum sources in the OMC-2/3 region was estimated to 
be several xlO^ 6-7 ^ cm -3 from the 3.3 mm continuum emission measured with the NMA 
(Takahashi et al. 2008a, Shimajiri et al. 2008). While the angular resolutions of most of 
the previous 3.3 mm continuum observations are somewhat lower than the SMA angular 
resolution, it is clear that the average hydrogen volume density toward MMS 6-main is 
roughly one order of magnitude higher than the value derived for MMS 7 and the other 
OMC-2/3 continuum sources. 

We also compared the MMS 6-main continuum source with the surrounding envelope 
properties. The average number density of the envelope gas detected in the H 13 CO + emission 
was estimated to be (1.6-5.0) xlO 7 cm -3 with a size scale of 5000 AU. This suggests that 
the detected compact continuum emission has approximately one order of magnitude smaller 
size and three order of magnitude higher volume- density as compared with the surrounding 
envelope structure detected in the H 13 CO + emission. The compact core detected in the 
interferometric continuum observations is brighter than surrounding parental core because 
of the high density region towards the center. 

Several of the intermediate- to high-mass protostellar core candidates are located in the 
OMC-1 and OMC-1 south (OMC-1S) regions. In those regions, several bright continuum 
sources have previously been detected with the high angular resolution of the SMA (Beuther 
et al. 2004, Zapata et al. 2005). Measured 0.9 mm total fluxes in the OMC-1 continuum 
sources (i.e., source I, source n, SMA 1, and Hot Core) by Beuther et al. (2004) have 
been reported to be 0.32 to 1.9 Jy. Measured 0.9 mm total fluxes of the strong continuum 
sources in the OMC-1S region range from approximately 0.3 to 0.8 Jy (L. A. Zapata, private 
communication based on Zapata et al. 2005). Hence, the comparison of the fluxes between 
MMS 6-main (F to tai=3.4 Jy at 0.9 mm) and those continuum sources in the OMC-1 and 
OMC-1S regions suggests that MMS 6-main is one of the brightest sources in the entire 
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OMC filament in terms of the total flux. 

The next interesting question is whether MMS 6-main is gravitationally bound or not. 
To discuss this issue, we compared the Virial mass with the mass derived from the inter- 
ferometric dust-continuum observations. The Virial mass is derived from the equation of 
M vlY ~5RCefi/G. Here, the effective sound speed of the gas is derived by C% s = AV r Q 2 bs /81n2, 
where G is gravitational constant, and R is radius of the envelope gas derived from the 
H 13 CO + observations. Using AV^bs^l-O km s _1 which was measured from the H 13 CO + (1-0) 
line profile at MMS 6-main peak with a beam size of 4". 6x3". 2 (1800x1300 AU), the virial 
mass was estimated to be 2.3 M . On the other hand, the H 2 mass derived from the 0.9 mm 
dust continuum emission is 3.0 M within a size scale of L>^510 AU. These results show 
that the H 2 mass derived from the 0.9 mm dust continuum emission is higher than the mass 
derived from the Virial mass estimate over a larger scale, suggesting that MMS 6-main is 
probably in the gravitationally bound state. 

In summary, MMS 6 is one of the brightest submilimeter continuum sources in the entire 
OMC region. The much higher hydrogen volume density, 6.5xl0 9 cm -3 , of MMS 6-main is 
consistent with its gravitationally bound nature. Its density is more than four to five orders 
of magnitude larger than the average hydrogen volume density of "typical star-forming dense 
cores", which is ^10 4-5 cm -3 . 



We probe the internal structure of the MMS 6-main core in order to improve our un- 
derstanding of the initial conditions of the star-forming core. In this section, we discuss the 
density distribution of the continuum source, MMS 6-main, using the visibility amplitude 
plot. The calibrated visibility data are given by the Fourier transform of the intensity dis- 
tribution of the observed target. In the Fourier domain analysis, a model of the core density 
structure can be directly compared with the interferometric data without the non-linear 
image deconvolution process and the missing flux problem. 

For optically thin dust emission, the observed intensity at impact parameter b is given 
by (e.g., Harvey et al. 2003a), 



where R out is the outer radius of the system, T^(r) is the dust temperature, p(r) is the 
density structure of the system, k v is the dust opacity (mass absorption coefficient), and 
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Internal Structure 
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B v is the Planck function. Here, we assume a spherically symmetric spatial structure. If 
the density and temperature follow simple radial power-laws, p(r)ocr~ p and T(r)ocr _g , and 
the opacity follows, k^ocz/^, the intensity has a simple power-law form, /(V, 6)oc^ 2+ ^6 _ ^ +g_1 ^. 
Since the attenuation resulting from the primary beam antenna pattern of the interferometer 
is negligible, the visibility distribution is also a power-law (i.e., Fourier transform of the 
intensity), V{y, ^)ocz/ 2+ ^^^ +g_3 ^. For reasons of expediency, hereafter, we define 7 = p+q— 3. 

Here, the temperature distribution is also important for extracting the density structure 
since the temperature and density profiles are coupled and observed as a single intensity 
profile. In the MMS 6-main case, the largest uncertainty comes from the unknown presence 
of a central protostar. We consider two kinds of temperature distribution: (i) isothermal 
temperature distribution: No heating source at MMS 6-main since no conclusive evidence 
for the central heating source of MMS 6-main was obtained from our observations and (ii) 
power-law temperature distribution: Assumption that the dust is optically thin in the 
envelope and the heating source (i.e., protostar) is located at the center. For the latter case 
where the optically thin dust envelope is heated by a central source, the dust temperature 
distribution is often expected to follow a power-law index determined by the radius and 
emissivity such as T^(r)ocL g / 2 r _g , where q = 2/(4 + /?), where L is the luminosity from the 
central star (Doty & Leung 1994). The emissivity value derived in Section 3.3, /3=0.93, for 
MMS 6-main was adopted, and q was estimated to be 0.41. 

Figure 7 shows the 3.3 mm and 0.9 mm binned visibility amplitudes as a function of the 
projected t^-distance. In order to discuss the density structure of MMS 6-main dusty core, 
two simple geometrical models, (i) Gaussian shaped density model and (ii) power-law density 
model, were compared with the observed visibility amplitudes with a baseline range up to 
80 kA (corresponds to the core radius of ^-450 AU). 

The 3.3 mm extended component is fitted fairly well by a 2D-Gaussian model with the 
size of 4 // .0x2 // .8 (1600x1100 AU; denoted by a solid curve in Figure 7a). On the other 
hand, the power-law model applied to the entire t^-range (denoted by a dashed curve) is 
definitely a poor fit. This can be seen to be due to the behaviors at the shortest baselines 
corresponding to the larger scale structures. At 3.3 mm, the emission does not rise sharply 
at the short t^-spacings as would be expected for a power-law distribution, consistent with 
a definite outer radius to the emission structures. A better power-law fit is derived for the 
3.3 mm visibility data restricted to larger baselines, >12 kA (see a doted curve). This can be 
explained by the absence of short spacing information to constrain the extended component. 
As in the case for the 3.3 mm data, the extended structure in the 0.9 mm data is fitted well 
with either the 2D-Gaussian model (with a size of 2 // .3x2 // .0, corresponds to 920x800 AU) 
or the power-law density distribution up to the uv-distance of 18 kA. These results suggest 
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that the extended halo around MMS 6-main has a definite outer radius of 2000-3000 AU 
(corresponds to 12-18 kA). 

The 3.3 mm extended component is fitted with a power-law index of 7=-0.55. Since 
p = 3 — q + 7, we estimate p=2A5 for the isothermal case (i.e, q=0) and p=2.00 for the 
power-law temperature case (i.e., g=0.41). The estimated density profile from the inter- 
ferometric data, with p=2.00~2.45, is roughly consistent with the large-scale (10000 AU) 
density structure, with p r ^2 ) as derived from previous single-dish 850 /xm observations by 
JCMT/SCUBA (J0rgensen et al. 2006). The observed density power-law index toward MMS 
6-main, therefore, roughly agrees with a quasi-static star-forming scenario with a singular 
isothermal sphere (SIS), which has a density distribution of p(r)ocr -2 (Larson 1969; Penston 
1969), or even has a steeper power-law index than the isothermal case. There is indication 
that a steeper density distribution is established as the spatial size-scale decreases toward 
450 AU. In the visibility fitting, the 0.9 mm data show a very similar amplitude power-law 
index (7=0.61) and hence similar density power-law index (p= 1.98-2. 39). 

One feature to note is that no flattened density structure was detected at the core 
center, which is usually expected toward prestellar cores @ (e.g., Ward-Thompson et al. 
1994). Non-detection of the flattened density structure implies two possibilities; (i) current 
angular resolution is still too coarse, as a centrally flattened structure might exist with a size 
scale much less than i?~450 AU (corresponds to 80 fcA), or (ii) MMS 6-main already has 
a central protostar, but the dynamically collapsing area, which has a density distribution 
of p(r)ocr -L5 such as proposed by Shu et al. (1977), is still compact (i.e., R <450 AU). In 
either case, a constant power-law index suggests that MMS 6-main is in either (i) the latest 
stage of prestellar core or (ii) the earliest stage of protostellar core. 

Observationally, the density profile maintains the same power-law index up to the uv- 
distance of 80 kA, which corresponds to the linear size of it/^-450 AU. This implies that the 
upper limit to the size of the inf ailing envelope is iiS-450 AU. C e ff^0.28 km s _1 as derived 
from the assumption of T^20 K, then implies an upper limit to the age of the protostar 
of r core <R/C e ff~7.6xl0 3 yr. On the other hand, the free-fall time of MMS 6-main as 
derived from the average density; r ff = l/ v / Gp-2.12xl0 5 (n H2 /10 5 cm- 3 )- 1 / 2 -830 yr. Both 
timescales are much shorter than the typical lifetime of the Class 0/1 protostars, which is 
^lO^ 4-5 ) yr (Beichman et al. 1986; Cohen, Emerson, & Beichman 1989; Kenyon et al. 1990; 
Beichman, Boulanger, & Moshir 1992). 



5 We have 7.3 mm data, which have iw-ranges from 4.5 to 477 kA. However, there are only sufficient 
sensitivities for baselines shorter than ~100 kA. Hence the 7.3 mm data did not add new information on the 
compact component, as suggested in the 0.9 mm and 3.3 mm data 
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In addition to the extended envelope, a compact source, which appears as a remnant 
unresolved component, was suggested from the Gaussian fitting of the 3.3 mm and 0.9 mm 
data (denoted by Figure 7a, b dashed line). Although the observed ^-coverage is not enough 
to clearly show this component, the constant or slowly decreasing visibility amplitudes at 
long baselines with -uu-distances of more than 60 kX (corresponds to a linear size scale of 
L>~1200 AU) suggest the presence of this compact source. In contrast, power law-density 
fittings do not require an unresolved component. Longer baseline data with a ^-distance 
more than 80 kA are necessary to judge whether the unresolved component exists or not. 
From the Gaussian fitting, the best-fitted offset flux is derived as 56 mJy at 3.3 mm and 1.6 
Jy at 0.9 mm, which suggests an upper limit of H 2 mass of 0.2-1.0 M Q under the assumption 
of Tj=20-100 K. 



5. Conclusion and Future Prospects: MMS 6-main is an Extremely Young 

Stage of the Protostellar Core? 

The millimeter and submillimeter multi-wavelength continuum observations with the 
CO (1-0) and H 13 CO + (1-0) line emissions have been performed with the NMA, SMA, and 
VLA. The main results are summarized as follows; 

• A compact continuum source, MMS 6-main was detected with an apparent H 2 mass of 
3.0 M Q and a typical size of 510 AU. Despite its compact appearance, no clear evidence 
of a radio jet or a CO molecular outflow were detected at the position of MMS 6-main. 
Moreover, from the mult i- wavelength continuum images, no infrared source is seen at 
a wavelength shorter than 8 fim at the peak position of MMS 6-main. The derived 
H 2 column density, 2.6 xlO 25 cm -2 , corresponds to the visual extinction of A^ 15000 
mag and A 8/im =790 mag. The estimated A 8/im toward MMS 6-main is much higher 
than those derived toward other YSOs in the OMC filament, ~10 mag. Therefore, it 
is possible that the central protostar associated with MMS 6-main cannot be detected 
because of the extremely high extinction. 

• At least two orders of magnitude higher volume density, ^6.5 xlO 9 cm -3 , as compared 
with the other OMC-2/3 continuum sources, was suggested toward MMS 6-main. The 
internal structure of MMS 6-main roughly agrees with an isothermal dusty core with 
a power-law index of p~2. We note that the density structure with p~2 continues 
up to the ^-distance of 80 kA, which corresponds to the linear size of iiS-450 AU. 
Assuming that the upper limit of the inf ailing radius is i?^450 AU, the upper limit 
of the protostellar age is estimated to be <7.6xl0 3 yr, while the free-fall time scale 
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of the core is estimated to be ^830 yr from the average dust density. Both timescales 
are roughly consistent with each other, and the estimated timescale is extremely short 
as compared with the typical time scale of the Class 0/1 protostars which is ^10^ 4-5 ^ 
yr. This suggests that MMS 6-main is probably in the pre-stellar phase or a transition 
phase from the pre-stellar to the proto-stellar evolutionary stage. An unresolved com- 
ponent was suggested from the Gaussian fitting of the 3.3 mm and 0.9 mm data with 
an upper limit to the mass of 0.2-1.0 M . The alternative power-law density fitting 
would not require such a compact component. Longer baseline data with a iw-distance 
more than 80 kA are necessary to constrain the inner most structure. 

Due to the lack of even higher angular resolution, sensitivity, and an appropriate molec- 
ular tracer, we could not directly measure the physical properties of the 100 AU scale cir- 
cumstellar disk or the possible central protostar. SMA and VLA observations with higher 
angular resolutions and also future ALMA observations will reveal the internal structure of 
MMS 6-main with a size-scale down to several AU. High angular resolution data of H 2 
masers could be useful for searching for the very beginning of the mass ejection phenomena 
(i.e., the micro-jet) from the central protostar such as for the S106 FIR case (Furuya et al. 
1999). Moreover, polarization observations with higher spatial resolutions can be useful for 
elucidating the magnetic field structure at this early evolutionary stage. In addition to the 
observational studies, comparison with models, which predict the very early stages of star- 
forming cores such as Bonner-Ebert spheres or adiabatic "First-Cores" (e.g., Larson 1969; 
Masunaga et al. 1998; Saigo 2006, 2008) will also be important. 
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Fig. 1. — (a); 850 /im continuum image taken with SCUBA/ JCMT (Johnstone et al. 1999). 
The contour intervals are 0.5 Jy beam -1 beginning at 0.5 Jy beam _1 . (b), (c), and (d); Map 
of 3.3, 2.3, and 7.3 mm continuum images in the MMS 6/5 regions (large scale). The contour 
levels start at 3a level, with an interval of 3a. The maps were made using all i^-coverage 
with robust weighting, (e), (f), (g), and (h); Focused map of 0.9, 3.3, 2.3, and 7.3 mm 

continuum images around MMS 6-main. The contour levels are -3, 3, 6, 10, 20, 30a, 

Negative contours are shown by dashed lines. The crosses show the peak positions of the 1.3 
mm continuum image identified by IRAM 30 m (Chini et al. 1997). Filled ellipses at the 
bottom right corners show the synthesized beams of each map. The detailed parameters of 
these maps are presented in Table 2. 
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Fig. 2. — (a) 2.2 //m-band and (b) 8.0 /im images overlaid with the 3.3 mm continuum maps 
made with the NMA (contours start at 3-a level with an interval of 6 a), (c) Focused image 
of Figure (b). Identified infrared sources are denoted by the dashed lines. Thin solid lines in 
the images show the boundary area for the infrared source identifications, which is roughly 
consistent with the distribution of the OMC filamentary structure. 
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Fig. 3. — Total intensity map of the CO (1-0) emission in contours, made with the NMA, 
superposed on the Xs-image in gray scale taken with SIRIUS/IRSF. The velocity ranges of 
the blueshifted and redshifted emission are Vlsr — 8.3 to 10.3 km s _1 and Vlsr = 12.5 to 
14.5 km s _1 , respectively. A filled ellipse at the bottom left corner shows the synthesized 
beam of the CO (1-0) emission. Open squares show peak positions of MMS 5, MMS 6-main 
and MMS 6-NE as detected by the interferometric data. Open circles show positions of 
identified infrared sources (i.e., IRS 1 to IRS 5) measured at Spitzer/IRAC bands. 
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Fig. 4. — (a) Total intensity map of the H 13 CO + (1-0) emission integrated from Vlsr — 
10.8 to 12.1 km s _1 (Large scale). The contour levels start at 2a level, with an interval of 
2a (l-a= 0.24 Jy beam -1 km s -1 ). Negative contours are shown by dashed lines. Open- 
squares and open-circles show positions of milllimeter/submillimeter sources; MMS 5/6; and 
infrared sources (i.e., IRS 1, 2, 3, 4, and 5), respectively. Filled ellipse at the bottom left 
corner shows the synthesized beam. Large circle on the map denotes the FWHM size of the 
primary beam. Figure (b) shows focused map of Figure (a). 
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Fig. 5. — Flux density of MMS 6-main as a function of the frequency. Black (filled and open) 
squares show total flux densities of each frequency measured from the images. The measured 
values are presented in Table 4. 8.4 GHz flux was referred to Reipurth et al. (1999; taken 
with VLA D-array observations) as an upper limit value. Each flux has an error bar within 
each square. The best fitting result is denoted by dashed line. 
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Fig. 6. — Dust emissivity index (/?) maps of the MMS 6 region (color) derived from the 3.3 
mm and 0.9 mm continuum maps overlaid with the 0.9 mm continuum image (contours). 
Both 3.3 mm and 0.9 mm continuum maps were convolved with the 3" beam size (see upper 
right corner), and values for (3 were calculated using the equation (1). 
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Fig. 7. — 4 kX binned visibility amplitudes as a function of ^-distance for (a) : 3.3 mm and 
(b) : 0.9 mm continuum data, together with the best-fit models. Gaussian+point sources 
fitting results are denoted by solid lines. Power-law fitting results are denoted by dashed 
lines. For Figure (a), two power-law fits; (i) the best power-law fit derived from all the 
visibility data was denoted by the dashed line, and (ii) the best power-law fit derived from 
the visibility data ranging from 12 kA to 80 kA (same as the SMA ^-coverage) was denoted 
by the dotted line. Horizontal dashed lines in the Figures show the best-fit offset values 
for the possible presence of an unresolved component, la statistical errors were denoted by 
asterisks with solid line. 



Table 1. Observational Parameters of Continuum Emissions 



Parameter 


u.y mm 


2.3 mm 


3.3 mm 


7.3 mm 


3.6 cm 


Telescope 


SMA 


NMA 


NMA 


VLA 


VLA 


Configurations 21 


compact 


D 


AB, C, and D 


C and DnC 


DnC 


Primary beam HPBW (arcsec) 


36 


46 


77 


60 


320 


Synthesized Beam HPBW (arcsec) 


2.2x2.0 (-32°) 


4.7x4.0 (16°) 


2.6x2.0 (-4.8°) 


1.4x0.6 (-87°) 


10x3.0 (67°) 


Projected base line range (kA) 


13-88 


4-38 


3 - 115 


4.5 - 477 


8.7 - 55 


Maximum detectable structure (arcsec) b 


13 


35 


56 


43 


180 


Bandwidth (GHz) 


4 


2 


2 


0.1 


0.1 


Gain calibrator 


0530+135 


0530+135 


0530+135 


05416-05418 


05416-05418 


Flux of the gain calibrator (Jy) 


1.7 


4.0 


2.5 


0.7 


1.4 


Bandpass calibrator 


3C273 and 3C454.3 


0530+135 


3C454.3 






Flux calibrator 


Titan 


Uranus 


Uranus 


3C286 and 3C48 


3C48 


System temperature in DSB (K) c 


200 - 400 


300 - 800 


240 - 600 


d 


d 


RMS noise level (mJy beam -1 ) 


43 


6.5 


1.0 


0.1 


0.05 - 0.06 



a D and AB are most compact and sparse configurations, respectively. 

b Our millimeter and submillimeter observations were insensitive to more extended emission than this size scale structure at the 10 
% level (Wilner Sz Welch 1994). The maximum detectable structure for the VLA observations were referred from the VLA home page 
( [http://www.vla.nrao.edu / astro / guides / vlas / current / nodel0.html ) . 

c System noise temperature of the SIS receiver was measured at the source elevation angle. 

d In stead of the system noise temperature information, VLA observation log including weather and system conditions is available on 
the VLA home page ( [http://www.vla.nrao.edu/operators/logs /2008/06/200 8-06-01-1758-AT359.pdfl > 
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Table 2. Observational Parameters of Line Emissions observed with the NMA 



Parameter 


H 13 CO+ (1-0) 


12 CO (1-0) 


Configuration* 1 


C and D 


D 


Baseline (kA) 


2.9 to 100 


7.8 to 30 


Primary beam HPBW (arcsec) 


77 


62 


Synthesized Beam HPBW (arcsec) 


4.6x3.2 (-23°) 


6.1x4.4 (-13°) 


Maximum detectable structure (arcsec) b 


60 


21 


Velocity Resolution (km s _1 ) 


0.22 


2.0 


Gain calibrator 


0531+135 


0531+135 


Bandpass calibrator d 


3C 454.3 


3C 273 


System temperature in DSB (K) e 


300 - 600 


200 - 600 


RMS noise level (Jy beam -1 ) 


9.2e-02 


2.0e-01 



a D configuration is the most compact configuration. 

b Our observations were insensitive to more extended emission than this size 
scale structure at the 10% level (Wilner & Welch 1994). 

c Gain calibrators were observed every 20 minutes. The flux density of the 
gain calibrator was determined relative to Uranus. 

d Bandpass calibrations were achieved by 30 to 40 minutes observations. 

e System noise temperature of the SIS receiver was measured at the source 
elevation angle. 



Table 3. Gaussian image fitting results toward MMS 6-main 



Wavelength 


R.A. (J2000) 


Dec. (J2000) 


Total flux 


Peak flux 


Deconvolved size a 


P.A. 


(mm) 






(mJy) 


(mJy beam -1 ) 


(arcsec) 


(degree) 


MMS 6-main 


0.9 mm 


5^35 m 23 s .42 


-5°01 / 30 // .54 


4.0e+03 


3.3e+03 


2.4x2.1 


-28 


2.3 mm 


5^35 m 23 s .42 


-5°01'30".61 


4.5e+02 


2.5e+02 


4.1x3.7 


172 


3.3 mm 


5^35 m 23 s .43 


-s^i'so^so 


1.4e+02 


8.0e+01 


2.2x1.8 


-31 


7.3 mm 


5^35 m 23 s .42 


-5°01 / 30 // .52 


6.0e+00 


4.1e+00 


0.8x0.5 


67 


MMS 6-NE 


0.9 mm b 


5^35 m 23 s .50 


-5°01 / 29 // .58 


1.6e+03 


3.2e+02 


5.5x3.8 


-15 


MMS 5 


3.3 mm 


5^35 m 22 s .44 


-5°01 / 14 // .50 


_b 


1.2e+01 


_b 


_b 


7.3 mm 


5^35 m 22 s .47 


-5 o 01'14".35 


b 


1.4e+00 


b 


b 



a l /; corresponds to 400 AU at the Orion Molecular Cloud. 

b Deconvolved size could not be derived from the 3.3 mm and 7.7 mm continuum data. 



Table 4. Identified infrared sources in the MMS 5/6 region 



Source Name 


RA. 


Dec. 


Note a 




(J2000) 


(J2000) 





IRS 1 (MMS 5) 5 35 22.4 -5 01 14.3 Class 

IRS 2 5 35 24.3 -5 01 20.8 PMS 

IRS 3 5 35 23.5 -5 01 29.0 Class I 

IRS 4 b 5 35 24.2 -5 01 31.2 Reflection nebula 

IRS 5 5 35 23.7 -5 01 40.0 Class I 



Evolutionary status of each source was determined from the 
Spectral Energy Distribution. 

b IRS 4 has extended emission, and its position was measured 
at the brightest part of the near-infrared structure. 



Table 5. Flux density of MMS 6 main 



Wavelength 


Flux density 


Error 


Deconvoled Size (P. A.) 


Om) 


(mJy) 


(mJy) 


(arcsec X arcsec) , (degree) 



3.6 cm a <1.5e-l 

7.3 mm 6.6e+00 b ' c 

3.3 mm 9.6e+01 b ' c 

2.3 mm 1.7e+02 b ' c 

0.9 mm 3.4e+03 



4.0e-02 — d 

1.9e-01 0.8x0.4 (160) 

3.5e+00 1.1x1.0 (106) 

2.6e+01 — d 

1.0e+02 1.0x0.8 (166) 



a Referred from Reipurth et al. (1999) 

b The flux densities were measured using uv-data which are longer than 
30 kA. 

c The values are different from the values derived from the gaussian 
image fitting (Table 3) because we applied the gaussian fitting using the 
different uv-range and also different method. 

d Deconvolved size could not be derived from the data. 



